Parton transport and hadronization from the dynamical quasiparticle point of view 
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The hadronization of an expanding partonic fireball is studied within the Parton-Hadron-Strings 
Dynamics (PHSD) approach which is based on a dynamical quasiparticle model (DQPM) matched 
to reproduce lattice QCD results in thermodynamic equilibrium. Apart from strong parton interac- 
tions the expansion and development of collective flow is driven by strong gradients in the parton 
mean- fields. An analysis of the elliptic flow 1)2 demonstrates a linear correlation with the spatial 
eccentricity e as in case of ideal hydrodynamics. The hadronization occurs by quark-antiquark fu- 
sion or 3 quark/3 antiquark recombination which is described by covariant transition rates. Since 
the dynamical quarks become very massive, the formed resonant 'pre-hadronic' color-dipole states 
(qq or qqq) are of high invariant mass, too, and sequentially decay to the groundstate meson and 
baryon octets increasing the total entropy. This solves the entropy problem in hadronization in a 
natural way. The resulting particle ratios turn out to be in line with those from a grandcanonical 
partition function at temperature T ~ 170 MeV rather independent from the initial temperature 
and indicate an approximate strangeness equilibration. 

PACS numbers: 25.75.-q, 13.60.Le, 14.40.Lb, 14.65.Dw 
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I. INTRODUCTION 

The 'Big Bang' scenario implies that in the first micro- 
seconds of the universe the entire state has emerged from 
a partonic system of quarks, antiquarks and gluons - a 
quark-gluon plasma (QGP) - to color neutral hadronic 
matter consisting of interacting hadronic states (and res- 
onances) in which the partonic degrees of freedom are 
confined. The nature of confinement and the dynamics of 
this phase transition has motivated a large community for 
several decades and is still an outstanding question of to- 
days physics. Early concepts of the QGP were guided by 
the idea of a weakly interacting system of partons which 
might be described by perturbative QCD (pQCD). How- 
ever, experimental observations at the Relativistic Heavy 
Ion Collider (RHIC) indicated that the new medium cre- 
ated in ultrarelativistic Au+Au collisions was interact- 
ing more strongly than hadronic matter (cf. [1] and Refs. 
therein) and consequently this concept had to be given 
up. Moreover, in line with theoretical studies in Refs. 
[2-4] the medium showed phenomena of an almost per- 
fect liquid of partons [5, 6] as extracted from the strong 
radial expansion and elliptic flow of hadrons [5] . 

The question about the properties of this (nonpertur- 
bative) QGP liquid is discussed controversially in the lit- 
erature and dynamical concepts describing the formation 
of color neutral hadrons from partons are scarce [7-13] . 
A fundamental issue for hadronization models is the con- 
servation of 4-momentum as well as the entropy problem 
because by fusion/coalescence of massless (or low con- 
stituent mass) partons to color neutral bound states of 
low invariant mass (e.g. pions) the number of degrees 
of freedom and thus the total entropy is reduced in the 
hadronization process [9-11]. This problem - a violation 
of the second law of thermodynamics as well as of the 
conservation 4-momentum and flavor currents - definitely 



needs a sound dynamical solution. 

A consistent dynamical approach - valid also for 
strongly interacting systems - can be formulated on the 
basis of Kadanoff-Baym (KB) equations [14, 15] or off- 
shell transport equations in phase-space representation, 
respectively [15-17]. In the KB theory the field quanta 
are described in terms of propagators with complex self- 
energies. Whereas the real part of the selfenergies can 
be related to mean- field potentials, the imaginary parts 
provide information about the lifetime and/or reaction 
rates of time-like 'particles' [4]. Once the proper (com- 
plex) selfenergies of the degrees of freedom are known the 
time evolution of the system is fully governed by off-shell 
transport equations (as described in Refs. [15-17]). 

The determination/extraction of complex selfenergies 
for the partonic degrees of freedom has been performed 
in Refs. [4, 18, 19] by fitting lattice QCD (1QCD) 'data' 
within the Dynamical QuasiParticle Model (DQPM). In 
fact, the DQPM allows for a simple and transparent in- 
terpretation of lattice QCD results for thermodynamic 
quantities as well as correlators and leads to effective 
strongly interacting partonic quasiparticles with broad 
spectral functions. We stress that mean-field potentials 
for the 'quarks' and 'gluons' as well as effective interac- 
tions have been extracted from 1QCD within the DQPM 
as well (cf. Ref. [19]). 



II. THE PHSD APPROACH 

The Parton-Hadron-String-Dynamics (PHSD) ap- 
proach is a microscopic covariant transport model that 
incorporates effective partonic as well as hadronic degrees 
of freedom and involves a dynamical description of the 
hadronization process from partonic to hadronic matter. 
Whereas the hadronic part is essentially equivalent to the 
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conventional Hadron-Strings-Dynamics (HSD) approach 
[20] the partonic dynamics is based on the Dynamical 
QuasiParticle Model (DQPM) [18, 19] which describes 
QCD properties in terms of single-particle Green's func- 
tions (in the sense of a two-particle irredicible (2PI) ap- 
proach) . 



With the spectral functions fixed by Eqs. (l)-(3) the 
total energy density in the DQPM (at vanishing quark 
chemical potential) can be evaluated as 

t °° =d »r^i^ 2 - 2 /^,pwc/t) (4) 



A. Reminder of the DQPM 

We briefly recall the basic assumptions of the DQPM: 
Following Ref. [21] the dynamical quasiparticle mass (for 
gluons and quarks) is assumed to be given by the thermal 
mass in the asymptotic high-momentum regime, which is 
proportional to the coupling g(T/T c ) and the tempera- 
ture T with a running coupling (squared), 
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Here N c — 3 stands for the number of colors while Nf 
denotes the number of flavors. The parameters control- 
ling the infrared enhancement of the coupling A = 2.42 
and T s = 0.46T C have been fitted in [21] to 1QCD results 
for the entropy density s(T). An almost perfect repro- 
duction of the energy density e(T) and the pressure P(T) 
from 1QCD is achieved as well. As demonstrated in Fig. 
1 of Ref. [19] this functional form for the strong coupling 
a s = g 2 /{Air) is in full accordance with the 1QCD calcu- 
lations of Ref. [22] for the long range part of the q — q 
potential, too. 

The width for gluons and quarks (for vanishing chem- 
ical potential p, q ) is adopted in the form [23] 

'■(?)■<'> 

where c = 14.4 (from Ref. [4]) is related to a magnetic 
cut-off. We stress that a non-vanishing width 7 is the 
central difference of the DQPM to conventional quasi- 
particle models [24-26]. It influence is essentially seen 
in correlation functions as e.g. in the stationary limit 
of the correlation function in the off-diagonal elements of 
the energy-momentum tensor T kl which defines the shear 
viscosity r\ of the medium [4, 27]. Here a sizable width is 
mandatory to obtain a small ratio in the shear viscosity 
to entropy density i]/s. 

In line with [21] the parton spectral functions thus are 
no longer S— functions in the invariant mass squared but 
taken as 
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separately for quarks and gluons (j = q,q,g). With the 
convention E 2 (p) = p 2 + M 2 - 7? , the parameters M 2 
and 7j arc directly related to the real and imaginary parts 
of the retarded self-energy, e.g. H, = M 2 — 2i-fjuj. 
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where tib and np denote the Bose and Fermi functions, 
respectively. The number of transverse gluonic degrees 
of freedom is d g — 16 while the fermic degrees of free- 
dom amount to d q = 4N c Nf = 36 in case of three flavors 
(Nf—3). The pressure P then may be obtained by inte- 
grating the differential thermodynamic relation 



p - 7 1 — = -T ao 
dT 

with the entropy density s given by 



s = 



dP _T ao + P 
df ~ f 
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This approach is thermodynamically consistent and rep- 
resents a 2PI approximation to hot QCD (once the free 
parameters in (1) and (2) are fitted to lattice QCD results 
as in Refs. [4, 18, 19]). 

As outlined in detail in Refs. [18, 19] the energy den- 
sity functional (4) can be separated in space-like and 
time-like sectors when the spectral functions aquire a fi- 
nite width. The space-like part of (4) defines a poten- 
tial energy density V p since the field quanta involved are 
virtuell and correspond to partons exchanged in interac- 
tion diagrams. The time-like part of (4) corresponds to 
effective field quanta which can be propagated within the 
light-cone. Related separations can be made for virtuell 
and time-like parton densities [18, 19]. Without repeat- 
ing the details we mention that mean-field potentials for 
partons can be defined by the derivative of the potential 
energy density V p with respect to the time-like parton 
densities and effective interactions by second derivatives 
of V p (cf. Section 3 in Ref. [19]). 



B. Hadronization in PHSD 

Based on the DQPM we have developed an off-shell 
transport approach denoted as PHSD where the degrees- 
of-freedom are dynamical quarks, antiquarks and gluons 
(q, q, g) with rather large masses and broad spectral func- 
tions in line with (1) - (3) as well as the conventional 
hadrons (described in the standard HSD approach [20]). 
On the partonic side the following elastic and inelastic 
interactions are included qq <-» qq, qq <-> qq, gg <-» gg, 
99 <7i 99 9 exploiting 'detailed-balance' with in- 
teraction rates from the DQPM [4, 18, 19] (cf. Section 
IIIB). The hadronisation, i.e. transition from partonic to 
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hadronic degrees of freedom, is described by local covari- 
ant transition rates e.g. for q + q fusion to a meson m of 
four-momentum p = (u, p) at space time point x = (t, x): 



as well as flavor current conservation, i.e. 
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X ^q Pq(Pq) ^q PqiPq) \Vqq\ 2 W m (x q ~ Xq,p q - p q ) 

xN q (x q ,p q ) Ng(xg,pq) £(flavor, color). (7) 
In (7) we have introduced the shorthand notation 

^=E / d 4 x j d 4 p j /(2n) 4 (8) 

where denotes a summation over discrete quantum 
numbers (spin, flavor, color); Nj(x,p) is the phase-space 
density of parton j at space-time position x and four- 
momentum p. In Eq. (7) ^(flavor, color) stands symbol- 
ically for the conservation of flavor quantum numbers as 
well as color neutrality of the formed hadron m which can 
be viewed as a color-dipole or 'pre-hadron'. Furthermore, 
Vqq(pp) is the effective quark-antiquark interaction from 
the DQPM defined by Eq. (31) and displayed in Fig. 10 
of Ref. [19] as a function of the local parton (q + q + g) 
density p p (or energy density). Furthermore, W m (x,p) is 
the phase-space distribution of the formed 'pre-hadron'. 
It is taken as a Gaussian in coordinate and momentum 
space (following Ref. [28]) with width V< r 2 > = 0.66 
fm (in the rest frame) which corresponds to an average 
rms radius of mesons. The width in momentum space 
is fixed by the uncertainty principle, i.e. AxAp = 1 (in 
natural units). We note that the final hadron formation 
rates are approximately independent on these parameters 
within reasonable variations. 

In principle, the two-particle Greens-function 
G<{Xq ,Pq,Xq,Pq) should appear in Eq. (7). The 
approximation of the two-particle Greens function by 
a (symmetrized/antisymmetrized) product of single- 
particle Greens functions is always a first step in a 
cluster expansion for Greens functions and neglects 
'residual correlations' stemming from higher order 
contractions. The same holds for an approximation 
of the three-particle Greens function by the (sym- 
metrized/antisymmetrized) product of single-particle 
Greens functions (cf. Ref. [29]). However, the DQPM 
with its dynamical spectral functions already includes 
the effects of strong two-body correlations - contrary to 
bare Green functions - such that the effect of residual 
interactions might be discarded in a first approximation. 
But there is no 'a priori' guarantee that this approxi- 
mation is appropriate under all circumstances. This, in 
principle, should be examined by lattice QCD in order 
to test the cluster decomposition in hot QCD. 

Related transition rates (to Eq. (7)) are defined for the 
fusion of three off-shell quarks (qi + q2 + q 3 <-> B) to a 
color neutral baryonic (B or B) resonances of finite width 
(or strings) fulfilling energy and momentum conservation 
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XUJ qi PqAPl) Uq 2 Pq 2 (p 2 ) ^<Z3 PqsiPs) 
x\M q qq\ 2 W B {x\ , X 2 , X 3 , Pi , p 2 , Ps) 

xN gi (xi,pi) N Q2 (x 2 ,P2) N g3 (x3,ps) ^(flavor, color) . 

Here the quantity Wb denotes the baryon phase-space 
distribution (evaluated in Jacobi coordinates) which is 
taken again of Gaussian shape with a width of 1 fm in 
coordinate space which corresponds to an average rms 
radius of excited baryons. The matrix element squared 
| Afgq g | 2 reflects the strength of three quark fusion pro- 
cesses and is fixed as follows: Since Regge trajectories 
for excited mesonic and baryonic states have the same 
slope (or string constant in the color dipole picture) we 



tentatively set \M q 



in our present work which 



implies that (so far) there is no need to introduce any new 
parameter. 

On the hadronic side PHSD includes explicitly the 
baryon octet and decouplet, the 0~ and 1~ meson 
nonets as well as selected higher resonances as in HSD 
[20]. Hadrons of higher masses (> 1.5 GcV in case of 
baryons and > 1.3 GeV in case of mesons) are treated 
as 'strings' (color-dipoles) that decay to the known (low 
mass) hadrons according to the JETSET algorithm [30]. 



III. HADRONIZATION OF AN EXPANDING 
PARTONIC FIREBALL 

We now turn to actual results from PHSD for the 
model case of an expanding partonic fireball at initial 
temperature T = 1.7 T c (T c = 0.185 GeV) with quasipar- 
ticle properties and four-momentum distributions deter- 
mined by the DQPM at temperature T = 1.7 T c . 



A. Initial conditions 

The initial distribution for quarks, antiquarks and glu- 
ons in coordinate space is taken as a Gaussian ellipsoid 
with a spatial eccentricity 



e=(y 2 -x 2 )/(y 2 +x 2 ) 



(10) 



and < z 2 >=< y 2 > in order to allow for the built-up 
of elliptic flow (as in semi-central nucleus-nucleus colli- 
sions at relativistic energies). In order to match the ini- 
tial off-equilibrium strange quark content in relativistic 
pp collisions the number of s (and s quarks) is assumed 
to be suppressed by a factor of 3 relative to the abun- 
dance of u and d quarks and antiquarks. In this way we 
will be able to investigate additionally the question of 
strangeness equilibration. 
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As mentioned above the dynamical evolution of the 
system is then entirely described by the transport dynam- 
ics in PHSD incorporating the off-shell propagation of the 
partonic quasiparticles according to Refs. [16] as well as 
the transition to resonant hadronic states (or 'strings') 
Eqs. (7), (9). The time integration for the testparti- 
cle equations of motion (cf. Refs. [16]) is performed in 
the same way as in case of hadronic off-shell transport 
where (in view of the momentum-independent width 7) 
the simple relation (19) in Ref. [31] is employed. For 
the collisions of partons two variants are at our disposal: 
i) geometrical collision criteria as employed in standard 
hadronic transport, ii) the in-cell method developed in 
Ref. [32]. The latter can easily be extended to describe 
2 <-> 3 processes etc. in a covariant way [33] and is the 
better choice at high particle densities (cf. Ref. [34]). 
The hadronization is performed by integrating the rate 
equations (7) and (9) in space and time which are dis- 
cretized by At and AV(t). We use local cells of volume 
dV(t)= 0.25 (1+bt) 3 fm 3 where t is given in units of fm/c 
and b= 0.7 c/fm. This choice approximately corresponds 
to a comoving grid for the expanding system. In each 
time-step At and cell AV the integrals in (7) and (9) are 
evaluated by a sum over all (time- like) testparticles using 
(e.g. for the quark density) 
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FIG. 1: (Color online) Upper part: Time evolution of the total 
energy E to t (upper line), the partonic contributions from the 
interaction energy V p and the energy of time-like partons T p 
in comparison to the energy contribution from formed mesons 
E m and baryons (+ antibaryons) E B+ g. Lower part: Time 
evolution in the parton, meson and baryon number for an 
expanding partonic fireball at initial temperature T = 1.7 T c 
with initial eccentricity e — 0. 



B. Dynamical evolution 



where the sum over J q implies a sum over all testparticles 
of type q (here quarks) in the local volume AV in each 
parallel run. In order to obtain lower fluctuations the 
integrals are averaged over the number of parallel runs 
(typically a few hundred). For each individual testpar- 
ticle (i.e. x q and p q fixed) the additional integrations in 
(7) and (9) give a probability for a hadronization process 
to happen; the actual event then is selected by Monte 
Carlo. Since energy-momentum conservation fixes the 
four-momentum p of the hadron produced - the space- 
time position x is fixed by (7) or (9) - the latter is rep- 
resented by a hadronic state with flavor content fixed 
by the fusing quarks (antiquarks). The latter decays to 
the lower mass hadrons according to JETSET [30] above 
thresholds of 1.3 GeV for mesonic states or 1.5 GeV for 
baryonic states (as in HSD). Lower mass hadrons (octet 
and decouplet states) are determined by the weight of 
their respective spectral functions at given invariant mass 
and selected by Monte Carlo. Note that the propaga- 
tion of partons includes the space-time derivatives of the 
quark and gluon mean-fields specified in Eq. (29) and 
displayed in Fig. 9 of Ref. [19]. 



In Fig. 1 (upper part) we show the energy balance for 
the expanding system at initial temperature T = l.7T c 
and eccentricity e =0, i.e. a fireball of spherical symme- 
try. The total energy E to t (upper line) - which at t = is 
given by (4) integrated over space - is conserved within 
3% throughout the partonic expansion and hadronization 
phase such that for t > 8 fm/c it is given essentially by 
the energy contribution from mesons and baryons (-(-an- 
tibaryons). The initial energy splits into the partonic 
interaction energy V p (cf. Eq. (19) in Ref. [19]) and the 
energy of the time-like (propagating) partons 

t p = J2 \/p? + M i(p P ) (i 2 ) 

i 

with fractions determined by the DQPM [19]. In Eq. 
(12) the summation over i runs over all testparticles in 
an individual run. The hadronization mainly proceeds 
during the time interval 1 fm/c < t < 7 fm/c (cf. the 
lower part of Fig. 1 where the time evolution of the q,q, g, 
meson and baryon (+antibaryon) number is displayed). 

As one observes from Fig. 1 on average the number 
of hadrons from the resonance or 'string' decays is larger 
than the initial number of fusing partons. This might be 
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FIG. 2: (Color online) The invariant mass distribution for fus- 
ing qq pairs (solid line) as well as qqq (and qqq) triples (dashed 
line) that lead to the formation of final hadronic states for an 
expanding partonic fireball at initial temperature T — 1.7 T c 
with initial eccentricity e = 0. 



astonishing since by partonic fusion the number of final 
states is conventionally reduced in coalescence models. 

In order to shed some light on the hadronization pro- 
cess in PHSD we display in Fig. 2 the invariant mass 
distribution of qq pairs (solid line) as well as qqq (and 
qqq) triples (dashed line) that lead to the formation of 
final hadronic states. In fact, the distribution for the 
formation of baryon (antibaryon) states starts above the 
nucleon mass and extends to high invariant mass cover- 
ing the nucleon resonance mass region as well as the high 
mass continuum (which is treated by the decay of strings 
within the JETSET model [30]). On the 'pre- mesonic 
side the invariant-mass distribution starts roughly above 
the two-pion mass and extends up to continuum states 
of high invariant mass (described again in terms of string 
excitations). The low mass sector is dominated by p, a±, 
uj or K*,K* transitions etc. The excited 'pre-hadronic' 
states decay to two or more 'pseudoscalar octet' mesons 
such that the number of final hadrons is larger than the 
initial number of fusing partons. 

Accordingly, the hadronization process in PHSD leads 
to an increase of the total entropy and not to a decrease 
as in case of coalescence models [9, 10]. This is a di- 
rect consequence of the finite (and rather large) dynam- 
ical quark and antiquark masses as well as mean-field 
potentials which - by energy conservation - lead to 'pre- 
hadron' masses well above those for the pseudo-scalar 
meson octett or the baryon octett, respectively. This 
solves the entropy problem in hadronization in a natural 
way and is in accordance with the second law of thermo- 
dynamics! 

The parton dynamics itself is governed by their propa- 
gation in the time-dependent mean-field U p {p p ) which is 
adopted in the parametrized form (as a function of the 
parton density p p ) given by Eq. (29) in Ref. [19]. Since 
the mean-field U p is repulsive the partons are accelerated 
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FIG. 3: (Color online) The effective gluon-gluon — » gluon- 
gluon (solid line), gluon-quark — > gluon-quark and quark- 
quark — » quark-quark (dashed line) cross section from the 
DQPM as a function of the parton density p p . 
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FIG. 4: (Color online) The interaction rate for the channels 
q + q — > g (dashed line, upper part) and g — > q + q (red 
line, upper part) for an expanding partonic fireball at initial 
temperature T = 1.7 T c with initial eccentricity e — 0. The 
lower part shows the collision rate of gluons (solid line) and 
quarks or antiquarks (dashed line). 



during the expansion phase on expense of the potential 
energy density V p which is given by the integral of U p 
over p p (cf. Section 3 in Ref. [19]). The interaction rates 
of the partons are determined by effective cross sections 
which for gg scattering have been determined in Ref. [4] 
as a function of T/T c . The latter are re-parametrized in 
the actual calculation as a function of the parton density 
using the available dependence of p P (T) on the tempera- 
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FIG. 5: (Color online) Upper part: The transverse momen- 
tum spectra 1/prdN /dpr for pions and nucleons for an ex- 
panding partonic fireball at initial temperature T — 1.7 T c 
with initial eccentricity e = 0. Lower part: The nucleon to 
pion ratio as a function of the transverse momentum pr cor- 
responding to the spectra in the upper part. The thin blue 
line is drawn to guide the eye due to the limited statistics of 
the PHSD calculation. 



ture T from the DQPM. 

The actual values for gg scattering are shown in Fig. 3 
as a function of the parton density p p (solid line) and 
demonstrate that gg cross sections up to 20 mb can 
be reached at p v « 2.5 fm~ 3 . The effective cross sec- 
tion drops rapidly with increasing p v which signals that 
weakly interacting partons might show up at very high 
parton density. Some note of caution has to be added 
here since though the cross section a gg drops with p p the 
collision rate of a gluon (~ cr gg p p ) increases slightly with 
Pp. For quark-quark or quark- antiquark elastic scatter- 
ing the cross section is reduced by a factor 4/9 in line 
with Casimir scaling. Quark-gluon elastic scattering (in 
the present implementation) is also reduced by a factor 
4/9 (cf. Fig. 3, dashed line). The channels qq — > g 
are described by a relativistic Breit-Wigner cross section 
which is determined by the actual masses of the fermions, 
the invariant energy y/s and the resonance parameters of 
the gluon (from the DQPM). In this case a further con- 
straint on flavor neutrality and open color is employed. 
The gluon decay to a uu, dd or ss pair is fixed by detailed 
balance. Further channels are gg <-*■ g which are given 
by Breit-Wigner cross sections (with the gluon resonance 
parameters) and detailed balance, respectively. 

The actual interaction rates for the channels q + q — > g 
and g — > q + q are displayed in Fig. 4 by the dashed line 



(upper part) and the solid line (upper part), respectively, 
for the expanding partonic fireball at initial temperature 
T = 1.7 T c with initial eccentricity e = 0. Within statis- 
tics the numerical result shows that detailed balance ac- 
tually is fulfilled for the expanding partonic system which 
was initialized in thermal equilibrium. The lower part of 
Fig. 4 shows the total number of collisions per time for 
gluons (solid line) and quarks or antiquarks (dashed line) 
which is higher for the fermions since the latter are much 
more frequent than the gluons. 

Let's, furthermore, have a look at the transverse mo- 
mentum pt spectra of hadrons emerging from the PHSD 
dynamics of the expanding fireball initialized with T = 
1.7 T c . The resulting px spectra are displayed in Fig. 5 
(upper part) for pions and nucleons and show that the 
nucleons become more abundant for px > 2.5 GeV/c. 
The ratio of nucleons to pions is depicted in the lower 
part of Fig. 5 and clearly demonstrate that baryons (an- 
tibaryons) become more frequent than mesons at high 
Pt > 2.5 GeV. This observation is in close analogy to the 
experimental findings in Au+Au collisions at top RHIC 
energies [5]. Note that a quantitative comparison with 
RHIC data is not meaningful due to the rather simpli- 
fied and special initial conditions employed here. 



C. Comparison to the statistical hadronization 
model 

It is, furthermore, interesting to have a look at the final 
particle ratios K + /ir + , p/tt + , A/K + etc. (after hadronic 
decays) which are shown in Table 1 . The latter ratios are 
compared to the grandcanonical statistical hadronization 
model (SM) [35-37] at baryon chemical potential p,B = 0. 
For p,B — the particle ratios depend only on temper- 
ature T and one may fix a freeze-out temperature, e.g., 
by the proton to 7r + ratio. A respective comparison is 
given also in Table 1 for T = 160 MeV and 170 MeV for 
the SM which demonstrates that the results from PHSD 
are close to those from the SM for T ps 170 MeV. This 
also holds roughly for the A/K + ratio. On the other 
hand the K + /tt + ratio only smoothly depends on the 
temperature T and measures the amount of strangeness 
equilibration. Recall that we initialized the system with 
a relative strangeness suppression factor of 1/3. The de- 
viation from the SM ratio by about 13% indicates that 
strangeness equilibration is not fully achieved in the cal- 
culations. This is expected since the partons in the sur- 
face of the fireball hadronize before chemical equilibra- 
tion may occur. A detailed discussion of results will be 
presented in a forthcoming study. 

The agreement between the PHSD and SM results for 
the baryon to meson ratio in the strangeness S=0 and 
S=l sector may be explained as follows: Since the final 
hadron formation dominantly proceeds via resonance and 
string formation and decay - which is approximately a 
microcanonical statistical process [38] - the average over 
many hadronization events with different energy/mass 





P/-K+ 


A/K+ 


K+/7V+ 


PHSD 


0.086 


0.28 


0.157 


SM T = 160 MeV 


0.073 


0.22 


0.179 


SM T = 170 MeV 


0.086 


0.26 


0.180 



TABLE I: Comparison of particle ratios from PHSD with the 
statistical model (SM) [37] for T= 160 MeV and 170 MeV. 



and particle number (in the initial and final state) leads 
to a grandcanonical ensemble. The latter (for hb — 0) is 
only characterized by the average energy or an associated 
Lagrange parameter (3 = l/T. 
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D. Elliptic flow 

Of additional interest are the collective properties of 
the partonic system during the early time evolution. In 
order to demonstrate the built up of elliptic flow we show 
in Fig. 6 the time evolution of 

V2 = {(pI-pI)/(pI+pI)) (13) 

for partons (solid line), mesons (long dashed line) and 
baryons (dot-dashed line) for an initial eccentricity e = 
0.33. As seen from Fig. 6 the partonic flow develops 
within 2 fm/c and the hadrons essentially pick up the 
collective flow from the accelerated partons. The hadron 
i>2 is smaller than the maximal parton V2 since by parton 
fusion the average vi reduces and a fraction of hadrons is 
formed early at the surface of the fireball without a strong 
acceleration before hadronization. We briefly mention 
that the reduction of the average hadron elliptic flow 
essentially is due to the finite parton masses which are 
larger than the temperature in the hadronization phase. 
This reduced hadron V2 is in contrast to the coalescence 
of massless partons. We discard a more detailed discus- 
sion and investigation of Eqs. (7) and (9) and refer the 
reader to an upcoming study. 

It is important to point out that in PHSD the elliptic 
flow of partons predominantly stems from the gradients 
of the repulsive parton mean-fields (from the DQPM) at 
high parton (energy) density. To demonstrate this state- 
ment we show in Fig. 6 the result of a simulation with- 
out elastic partonic rescattering processes by the short 
dashed line. 

Fig. 7 shows the final hadron V2 versus the initial ec- 
centricity e and indicates that the ratio V2 /e is practically 
constant (w 0.2) as in ideal hydrodynamics (cf. Fig. 3 
in Ref. [39]). Accordingly the parton dynamics in PHSD 
are close to ideal hydrodynamics. This result is expected 
since the ratio of the shear viscosity r/ to the entropy 
density s in the DQPM is on the level of r\js sa 0.2 [4] 
and thus rather close to the lower bound oirj/s = l/(47r) 
[40] . Note that the ratio 77/ s is dominantly determined by 
the quasiparticle width 7 (2) and low ratios on the level 
of rj/s « 0.2 require broad parton spectral functions as 
employed in the DQPM. 



FIG. 6: (Color online) Time evolution of the elliptic flow V2 
for partons and hadrons for the initial spatial eccentricity e = 
0.33 for an expanding partonic fireball at initial temperature 

r = i.7T c . 




e 

FIG. 7: (Color online) The hadron elliptic flow V2 versus the 
initial spatial eccentricity e for an expanding partonic fireball 
at initial temperature T = 1.7 T c . 

A further test of the PHSD hadronization approach is 
provided by the 'constituent quark number scaling' of the 
elliptic flow V2 which has been observed experimentally in 
central Au+Au collisions at RHIC [5, 41]. In this respect 
we plot V2/n q versus the transverse kinetic energy per 
constituent parton, 



with rriT and m denoting the transverse mass and actual 
mass, respectively. For mesons we have n q — 2 and for 
baryons/ antibaryons n q = 3. The results for the scaled 
elliptic flow are shown in Fig. 8 for mesons (dashed line) 
and baryons (solid line) and suggest an approximate scal- 
ing. We note that the scaled hadron elliptic flow V2/n q 
does not reflect the parton V2 at hadronization and is sig- 
nificantly smaller. Due to the limited statistics especially 
in the baryonic sector with increasing px this issue will 
have to be re-addressed with high statistics in the actual 
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FIG. 8: (Color online) The elliptic flow «2 - scaled by the num- 
ber of constituent quarks n q - versus the transverse kinetic 
energy divided by n q for mesons (dashed line) and baryons 
(solid line). 



heavy-ion case where the very early parton pt distribu- 
tion also shows 'power-law' tails. 

IV. SUMMARY 

In summary, the expansion dynamics of an anisotropic 
partonic fireball is studied within the PHSD approach 
which includes dynamical local transition rates from par- 
tons to hadrons (7), (9) (and vice versa). It shows col- 
lective features as expected from ideal hydrodynamics in 
case of strongly interacting systems. The hadronization 
process conserves four-momentum and all flavor currents 
and slightly increases the total entropy (by about 15% 
in the model case investigated here) since the 'fusion' 
of rather massive partons dominantly leads to the for- 
mation of color neutral strings or resonances that decay 
microcanonically to lower mass hadrons. This solves the 
entropy problem associated with the simple coalescence 
model! 

We find that the hadron abundancies and baryon to 
meson ratios are compatible with those from the sta- 



tistical hadronization model [35, 36] - which describes 
well particle ratios from AGS to RHIC energies - at a 
freeze-out temperature of about 170 MeV. Furthermore, 
strangeness equilibration is approximately achieved in 
the dynamical expansion and driven by the processes 
QQ 9 ss, which is a resonant process in the DQPM. 
However, although the final hadron ratios are compati- 
ble with a fixed freeze-out temperature (~ 170 MeV) we 
observe that the actual hadronization occurs at very dif- 
ferent energy densities (or temperatures) (cf. also Ref. 
[42]) such that the microscopic studies do not support 
the sudden freeze-out picture. 

Our calculations show that the hadron elliptic flow is 
essentially produced in the early partonic stage where 
also the strong repulsive parton mean-fields contribute 
to a large extent. This might explain why the hadron 
V2 from HSD calculations [43] underestimated the RHIC 
data on v-i essentially at midrapidity (in the pure 
hadron/string approach). The hadron elliptic flow from 
PHSD is smaller than the parton vi due to a partial can- 
cellation of the v 2 values from the individual partons in 
the fusion process to hadrons (but larger than in HSD). 
This cancellation essentially is a consequence of parton 
masses that are larger than the 'local temperature' during 
the hadronization phase. Nevertheless, our observations 
indicate an approximate quark number scaling (cf. Fig. 
8) for low and moderate transverse kinetic energies within 
the statistics reached so far. This issue will be followed 
up in more detail in a forthcoming investigation. 

The present study, however, serves only as a model case 
which allows for a more transparent interpretation of the 
various results. An application of the PHSD approach to 
ultrarelativistic heavy-ion collisions especially in compar- 
ison to differential experimental data is expected to shed 
further light on the transport properties of the partonic 
phase and the dynamics of hadronization. 

The authors like to thank A. Andronic for providing 
the SM results in Table 1 and for a critical reading of the 
manuscript. Furthermore, they are grateful to O. Linnyk 
for valuable discussions. 
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